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Abstract. We propose here a new approach for finding harmonic functions of killed 
random walks with small steps in the quarter plane. It is based on solving a certain 
functional equation that satisfies the generating function of the values taken by the 
harmonic functions. As a first application of our results, we obtain a simple expression 
for the harmonic function that governs the asymptotic tail distribution of the first exit 
time from the quarter plane. As another corollary, we prove, in the zero drift case, the 
uniqueness of the harmonic function. 



1. Introduction 

Context. Random processes conditioned on staying in cones of Z rf arouse a great interest 
in the mathematical community, as they appear in several distinct domains: quantum 
random walks [1, 2], random matrices [8], non-colliding random walks [5, 6, 9], etc. A usual 
way to realize this conditioning consists in using Doob /i-transforms, thanks to functions 
which are harmonic for the process, positive within the cone and equal to zero elsewhere — 
or equivalently harmonic and positive for the underlying killed process. It is therefore 
natural to be interested in finding some (all, if possible) positive harmonic functions for 
processes in cones of Z rf killed at the first exit time from the cone. 

In the literature, results on harmonic functions for killed processes in cones are broken 
up, and most of them concern very particular cones, as half spaces Z + x Z rf_1 and orthants 
Z^_. Regarding random walks with non-zero drift, general results are obtained in [16, 18] 
for the above domains when d = 2: the Martin boundary is found, and happens to be 
composed of infinitely many harmonic functions. For random walks with zero drift, the 
results are rare, and typically require a strong underlying structure: the random walks are 
cartesian products in [20]; they are associated with Lie algebras in [1, 2]; certain reflexion 
groups are supposed to be finite in [22]. This is problematic, as these harmonic functions 
are useful to construct other important processes (see above). Last but not least, knowing 
the harmonic functions for zero drift random walks in Z'Jr 1 is necessary for building the 
harmonic functions of walks with drift in Z£ , see [14]. 

In this article, we prove that for the whole class of walks with small steps and zero drift 
killed at the boundary of the quadrant Z+, there exists exactly one harmonic function. 
Further, the unique harmonic function is expressed in terms of sine and arcsine functions. 
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Our approach is new — to the best of our knowledge — and is based on functional equations 
that satisfy the generating functions of the values taken by the harmonic functions. 

Before presenting it (in point 5 below), we notice that in the literature, there are several 
ways to obtain harmonic functions: 

1. The most elementary approach just uses the definition of a harmonic function: if P 
denotes the generator of the Markov chain under consideration, we must have Pf = f 
inside of the cone, and / = outside. Unfortunately, except for a few particular cases (see, 
e.g., [26]), solving directly this discrete Dirichlet problem is out of range. 

2. A second one, based on Martin boundary theory, is used in [1, 14, 15, 16, 18, 22]. 
Briefly said, given a (transient) Markov process, this theory provides explicit expressions (as 
integral formulations) for all its (super)harmonic functions. To find the Martin boundary, 
it suffices [7] to compute the asymptotic behavior of quotients of Green functions. More 
details are given in Section 5.1. The main and profound difficulty of this approach lies in 
finding the (asymptotic of quotients of) Green functions. 

3. A third one, which is studied in [5, 6, 9, 17], is based on the fact that in many examples, 
the asymptotic tail distribution of the exit time t x from the cone of the random walk started 
at x is governed by a certain positive harmonic function V{x): as n — > oo, 



where x > 0. In comparison with point 2, this way presents the drawback of giving at 
most one harmonic function. 

4. Another possibility consists in comparing the solutions to the discrete and the classical 
(continuous) Dirichlet problems, and to use the numerous results which exist concerning 
the latter. Unfortunately, this comparison only exists for bounded domains (among which 
truncated cones), see [25]. 

5. A last and new approach consists in showing that the generating function of the values 
taken by the harmonic functions satisfies some functional equation, and in solving the 
latter. While the functional equation is easily obtained (it turns out to be a rewriting of 
the Dirichlet problem of point 1), its resolution will be more arduous! 

We now briefly present approach 5 in the one-dimensional case. Let X = (X(n)) ne z + 
be a random walk on Z + , with jumps to the two nearest neighbors, and killed 

at the boundary of Z + (the boundary is here simply the origin {0}). Let / be a harmonic 
function for this process. We have /(0) = and, for any iq ^ 1, 



Let H(x) = J2i >i /(*o)^ 0_1 - Multiplying (1.2) by x io ~ l and summing w.r.t. iq 1 yields 

(x 2 p_i - x + pi)H(x) = pi/(l). 



We wish to thank D. Denisov and V. Wachtel for suggesting us this new approach, during a workshop 
at Eindhoven University in February 2012. 



(1.1) 




(1.2) 
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It is then obvious to deduce an expression of H(x), and also, for any iq 1, of f(io), via 
Cauchy's formulae. 

Results. This work aims at presenting the approach summarized in point 5 (i.e., finding 
harmonic functions by solving functional equations) for the walks with small steps (i.e., 
with jumps only to the eight nearest neighbors) in the quarter plane Z?., see Figure 1. 
Our motivations to restrict ourselves to this class are the following: First, it is of special 




FIGURE 1. Walks with small steps in the quarter plane 

interest in probability [1, 4, 10, 11, 12, 16, 18, 22] as in combinatorics [3, 21]. Second, 
our way to solve the functional equation we shall obtain for (the generating function of) 
the harmonic functions needs complex analysis, and thus dimension 2. The hypothesis 
on the small steps in less crucial, and will be commented in Section 6. Further, we shall 
mostly assume that the increments of the random walk have zero mean, as it is the most 
interesting case (see above). This article is hinged on the following main results: 

1. We first state a functional equation satisfied by the generating function of the harmonic 
functions (Section 2). We solve it (Section 4), and we show that the solutions are closely 
related to a certain conformal mapping (introduced in Section 3). 

2. We then prove that for the class of zero- mean random walks with small steps in the 
quarter plane, there is exactly one harmonic function (Section 5). 

3. We also obtain the following two important corollaries (Section 5): in [6], it is proved 
that for the zero-mean random walks in the quarter plane, Equation (1.1) holds, where V(x) 
is shown to be harmonic, but is not made explicit. Here we find an expression for V(x), as 
a Taylor coefficient of a (relatively) simple function involving sine and arcsine functions. 
The second consequence is about lattice path enumeration, and more precisely about the 
counting of walks with small steps in the quarter plane. Recently, in [6, Section 1.5], 
the asymptotics of the number of excursions (walks starting and ending at given points) 
was obtained, up to some multiplicative constant involving two harmonic functions. The 
analysis we lead here allows us to obtain explicit expressions for both these functions. 

4. Finally, we present possible extensions of our approach and results (Section 6). Firstly, 
we show that our methods do work for random walks with non-zero mean increments. In 
particular, we obtain expressions for the harmonic functions in that case (this improves 
some results in [16, 18]). Regarding the non-zero drift case, we also prove that all harmonic 
functions converge to the unique harmonic function, as the drift goes to 0. Further, we see 
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that our approach also gives expressions for i-harmonic functions (i.e., functions / such 
that Pf = tf within the cone, P being the generator of the random walk and t > 0). 

2. Functional equations for harmonic functions 

2.1. Harmonic functions for walks with small steps in the quarter plane. Denote 
by (X,Y) = {X(n),Y(n)) n EZ + a random walk in the quarter plane Z^_, and let Pu j Q \ \S ] 
be the probability of event $ conditional on (X(0),Y(0)) = (io,jo)- Throughout we shall 
make the following assumptions: 

(HI) The walk is homogeneous inside of the quarter plane, with transition probabilities 
{pi,j}-i^i,j^.i to the eight nearest neighbors (po,o = 0); 

(H2) In the list pi : i,pi : o,pi,-i,Po,-ljP-l,-ijP-i,Ojf?-i,l>Po,lj there are no three consec- 
utive zeros; 

(H3) The drifts are zero: ipij = and jpi,j = 0. 

With assumption (HI) we can use the general framework for walks in the quarter plane 
developed by Fayolle, Iasnogorodski and Malyshev [11]. 

In this note we are interested in functions / which are harmonic for these random walks, 
i.e., in functions / = {f(io,jo))( io j ) e z 2 + sucn tliat: 

(PI) For any i ,j ^ 1, f(i ,h) = Y,-mi,j^iPi,jf( i + hh+j)- 
More specifically, we restrict ourselves to the harmonic functions which are 

(P2) zero on the boundary: if IqJq = then f(io,jo) = 0; 

(P3) positive inside of the quadrant: if io,jo > then /(io, Jo) > 0. 

A harmonic function / satisfying (P2) and (P3) is positive harmonic for the random walk 
killed at the boundary of the quarter plane. 

2.2. Structure of the remaining of Section 2. We fix a harmonic function / satisfying 
(PI), (P2) and (P3). We first obtain, in Section 2.3, a functional equation for the generating 
function 

(2-1) H(x,y) = E^i/^oJoK^V - 1 - 

This functional equation is stated in (2.4), and involves a certain bivariate polynomial 
L(x,y), see (2.3), that we study in Section 2.4. We then prove, in Section 2.5, that both 
generating functions 

(2.2) H(x,0) = E^i/^o,!)^ - 1 , H(0,y) = £ io3 * /(l, jo)^ " 1 , 

satisfy certain (rather simple) boundary value problems. Finally, starting from the latter, 
we obtain explicit expressions of H(x, 0) and H(0, y), then of H(x, y) via (2.4), and finally 
of f(io,jo), for any io, jo ^ 1) via t ne classical Cauchy's formulas. 

2.3. A first functional equation. The following polynomial — also called the kernel of 
the random walk — will be of the highest importance: 

(2-3) L(x,y) = xylE^ij^PijX-iy-i - 1]. 
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Lemma 1. For any random walk with property (HI) (but not necessarily (H2) and (H3) ), 
the generating function H(x,y) defined in (2.1) satisfies the functional equation 

(2.4) L(x, y)H(x, y) = L(x, 0)H(x, 0) + L(0, y)H(0, y) - L(0, 0)#(0, 0). 

Proof. The proof of (2.4) is a simple consequence of (PI) and (P2). □ 

We now make a series of remarks on this functional equation. 

Example. Let us verify (2.4) for the simple random walk, i.e., for the model with transition 
probabilities po,i = Pi = Po~i — P-l = 1/4. In that case, it is well known (see, e.g., [20]) 
that there exists a unique positive harmonic function (up to the multiplicative constants), 
with the product form /(in, Jo) = Wo- With (2.1) and (2.2), we then obtain 

(2-5) ^,0) = ^, H^y) = ^ w , H(x,y)= {i l _ x)2{i l _ y)2 . 

Furthermore, L(x,y) = (y/4)(x - l) 2 + (x/4)(y - l) 2 , L(x,0) = x/4, L(0,y) = y/A and 
L(0,0) = 0, see (2.3), and it becomes straightforward to verify that (2.4) does hold. 

Around the functional equation (2.4). In the literature, there actually exist many papers on 
problems that can be reduced to solving functional equations closed to (2.4). This typically 
arises in probability theory [10, 11, 22], in queueing theory [4], and in combinatorics as well 
[3, 21]. In our case, contrary to most of the functional equations present in the above-cited 
works: 

• The coefficient L(x, 0) (resp. L(0,y)) in front of unknown H(x,0) (resp. H(0,y)) 
depends only on one of the two variables x and y; 

• In the right-hand side of (2.4), there are no non-homogeneous terms depending 
on x and y. 

For these two reasons, the analytic description of the solutions to the functional equation 
(2.4) we shall propose here will be more precise than in other works. 

Domain of validity of (2.4). Equation (2.4) holds (at least) on the domain {(x,y) € C 2 : 
\x\ ^ 1, \y\ ^ 1}. Indeed, we have: 

Lemma 2. For a random walk satisfying to (HI), (H2) and (H3), the radius of convergence 
of H(x,0) and H(0,y) is one. Further, for any non-zero harmonic function, H(1,0) = 
#(0,1) = oo. 

The proof of Lemma 2 we shall give here 2 uses large deviation theory. It is worth noting 
that it is the only place in this article where a result is not derived as a consequence of the 
functional equation (2.4) (and it is an open problem to decide whether it is possible to do 
so). It is postponed to Appendix A. 

It is worth remarking that for a random walk (satisfying to (HI) and (H2) but) having a 
non-zero drift, Lemma 2 does not hold anymore. This is precisely this phenomenon which 
eventually will entail that the Martin boundary of a random walk with a non-zero drift is 
composed of infinitely many harmonic functions [16, 18]; see Section 6.1 for more details. 



'We are grateful to I. Ignatiuk-Robert for suggesting us the ideas of this proof. 
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2.4. Basic properties of the kernel. This part is purely technical, and presents some 
properties of the kernel L(x,y) introduced in (2.3), that are crucial for the sequel. 

Branch points. The kernel L(x,y) can also be written as 3 

(2.6) L(x, y) = a(x)y 2 + /3(x)y + -y(x) = a(y)x 2 + (3(y)x + j(y), 

where 
(2.7) 

a(x) =p_i _i£ 2 +p ,-i£+Pi P{x) =p_i i0 x 2 -x+pi i0 , t(x) =p„i i ix 2 +p ,i a; +Pi,i> 

a(y) =p-i _iy 2 +p_i, y+P-i,i, /3(y) =po ,-iy 2 - y+Po,i, i{y) = pi-iy 2 +pi,oy+pi,i- 

We also define 

(2.8) 5(x) = p(x f - 4a(x) 7 (x), %) = ^(y) 2 - 45(y)7(y). 

The roots of these polynomials are called the branch points of the kernel. The following 
facts regarding the polynomial 5 are proved in [11, Chapter 2]: Under (HI), (H2) and (H3), 
5 has degree three or four. Further, 1 is a root of order two of 5. We denote its roots by 

{^}i^ai witn 

(2.9) \xi\ ^ 1 ^ \x^\, X2 = X3 = 1, 

and X4 = 00 if 5 has degree three. We have x\ G [—1, 1) and X4 G (l,oo)U{oo}U(— 00, — 1]. 
Further, on R, 5(x) is negative if and only if x G [xi,X4] \ {1}. The polynomial 5 in (2.8) 
and its roots {ye}i^i^ satisfy similar properties. 

Algebraic functions defined by the kernel. In what follows, we denote by X(y) and Y(x) 
the algebraic functions defined by L(X(y),y) = and L(x,Y(x)) = 0. With (2.6) and 
(2.8) we have 

Functions X(y) and Y(x) both have two branches, called Xq, X\ and Yq, Y\, which are 
meromorphic on C \ [xi, X4] and C \ [2/1,2/4], respectively. We fix the notations by requiring 
that on the whole of C, 

\X (y)\ ^ |JTi(i/)|, \Y (x)\ < |Yi(x)|. 

On a curve defined by the complex values of function X(y). To conclude Section 2.4, we 
introduce the curve 

X([y u l]) = X ([ yi ,l])U Xt([ yi ,l]), 

with is symmetrical w.r.t. the real axis (since for y G [yi, 1], -^o(y) an d X\(y) are complex 
conjugates), and goes around point x\ (see [11, Theorem 5.3.3]). 



3 For readers used to [11], we notice that the notations we take here are not the usual ones, as the kernel 
L(x, y) in (2.3) is the "reversed" one of the kernel Q(x, y) in [11]: L(x, y) — x 2 y 2 Q(l/x, 1/y). 
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The curve X([yi, 1]) is smooth except at = 1, where it may have a corner point. 

More specifically, at 1 the angle between the segment [xi, 1] and the curve is given by [12] 



Denote by &x the open set bounded by the curve X([y\, 1]), which in addition contains 
the point x\ (or equivalently 0). For instance, in the case of the simple random walk (with 
P0,l = Pi,0 = P0,-l = P-1,0 = 1/4), X([y\, 1]) is the unit circle (see [11, Theorem 5.3.3]), 
hence is the unit disc. In that case, X([yi, 1]) is smooth at 1, which is coherent with 
9 = arccos(O) = vr/2 in (2.11). 

2.5. A second functional equation. Let <£x be the set introduced above, and let &x 
be its closure (in C U {oo}). This part aims at showing that H(x, 0) satisfies the following 
boundary value problem (below and throughout, x stands for the conjugate value of x £ C): 

(i) H(x,0) is analytic in <£x\ 

(ii) H(x,0) is continuous on ^x \ {1}; 

(iii) For all x on the boundary of ^x except at 1, L{x, 0)H(x, 0) — L(x, 0)H(x, 0) = 0. 

The problem of finding functions satisfying (i), (ii) and (iii) is a particular instance of a 
boundary value problem with shift (the complex conjugation plays in (iii) the role of the 
shift), see [19] for an extensive treatment of this topic. 

Proof of (i). Note that from Lemma 2, we already know that H(x, 0) is analytic within the 
open unit disc. Accordingly, if the domain &x is included in the unit disc, there is nothing 
to prove. Note that like the simple random walk (see Section 2.4), many walks do satisfy 
this property. 

On the other hand, there also exist walks such that part of ^x lies outside the unit disc. 
This is for example the case for the (Gessel's) walk with p\^\ = = P-i-i = P-i,o = 1/4, 
see [21]. In that case, (i) follows from an analytic continuation argument based on (2.4). 
This argument is exactly the same as that used in [21, Section 4], and we refer to this work 
for the details. □ 

Proof of (ii). See [21, Section 4], where a stronger statement (on the analytic continuation 
of the whole of C \ [1, X4]) is proved for a quite similar functional equation. □ 

Proof of (iii). Proving (iii) is classical [11, 18, 21]: it suffices to evaluate (2.4) both at Xo(y) 
and X\(y), for all y £ [y\, 1]. This way, the kernel L(x,y) vanishes, and the left-hand side 
of (2.4) too. Finally, taking the difference of the equations corresponding to Xo(y) and 



2.6. A key lemma. We conclude Section 2 by stating a lemma (all-present in Section 4) 
on a certain boundary value problem closed to that of Section 2.5. 

Lemma 3. Let h(x) be a function satisfying (i), (ii) and (iii), where 

(ii) h(x) is continuous on£?x- 
Then h(x) necessarily is a constant function. 



(2.11) 




Xi(y) yields (iii). 



□ 
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Proof. The proof of this result is contained in both [11] and [19]. 

3. A PERFECT CONFORMAL MAPPING 



□ 



In this section we introduce a certain function w(x) that will be crucial for our analysis. 
Indeed, we shall see that any function H(x, 0) solution to (i), (ii) and (iii) can be expressed 
as S(w(x))/L(x, 0), where S is a rational function and L(x,0) has been defined in (2.3) . 
Unexpectedly, this function w(x) happens to be a conformal mapping of a certain domain 
naturally attached to the random walk, namely, &x (introduced in Section 2.4). Luckily, 
this function w(x) has been carefully analyzed in [12]; in particular, an expression for w(x) 
was obtained. Let 



(3.1) /(*) 

and 
(3-2) 



d"(x A ) d'(x 4 ) 



6 

d"(0) 



6 



+ 



X — X4 

d"'(0)x 
6 



if Xi 7^ oo 
if X4 



T(x) 



oo, 



2/0*0 
d"(l) 



-1/2 



w(x) = sin ( — arcsin(T(x)) — — 
V u L 2 



with 6 as in (2.11). It is proved in [12] that w(x) enjoys the following nice properties: 
(Wl) w is analytic in Sfxi 
(W2) w is continuous on Sfy \ {!}; 

(W3) w establishes a conformal mapping of onto the complex plane C cut along a 
segment; 

(W4) There exists c / such that for x in the neighborhood of 1, 

c + o(l) 



(3.3) 



w(x) 



(l-ar)*/*' 

(W5) For all x on the boundary of Sfy, w(x) is real and w(x) = w(x). 
A priori, the function w(x) in (3.2) is well defined only for \T\ = \T(x)\ ^ 1. For |T| ^ 1, 
we have arcsin(T) = 7r/2±iln(T + VT 2 — 1), and as in [12, Section 2.2] we can write 



(3.4) 



w(x) 



T + \JT 2 - 1 



2tt/ 



2+ [T- \JT 2 - 1 



27r/0' 



Remark. It is also shown in [12] that w(x) is the unique function for which (W1)-(W5) 
hold (up to additive and multiplicative constants). 

Example. For the simple random walk, it is easy to compute X4 = 3+2^/2, d'(x±) = 4+3\/2, 
d"(x A ) = 11/2 + 3V2, d"(l) = -1/2, T(x) = [2(2 + y/2)(x- l)/(x-3-2 v / 2)]" 1/2 . Further, 
tt/9 = 2, so that with (3.4), 

x - 3 + 2\/2)(a; - 3 - 2^2) _ 1 2x 

2(x- l) 2 ~ 2 ~ (x-1) 2 ' 



(3.5) u;(x) = -4T 2 (T 2 - 1) 



4 We would like to insist on the fact that the function S(w(x))/L(x, 0) will not be the generating function 
associated with a harmonic function for any rational function S; however, it will solve the problem (i), (ii) 
and (iii) for any S. 
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4. Solving the functional equations 

In this part we characterize the harmonic functions satisfying (PI), (P2) and (P3). 
Precisely, we shall prove the following result: let w(x) be the function defined in Section 3. 
Subtracting by w(0), we assume that w(0) = 0. Define 

po,i/(M) 



w'(0) 



if pi t x = 0, 
if Pi,i/0. 



(4.1) u = -w(X (Q)), /i 

" w(X (0)) 

Theorem 4. Define \x and v as in (4.1). We have 

(4,, ^.»)=^- 

Of course, a similar expression could be written for H(0,y), after a suitable change of 
the parameters (namely, pij — > Pji). As for H(0, 0), it is simply equal to /(l, 1), see (2.1), 
and an expression for the complete generating function H(x,y) then follows from (2.4). 
Finally, the Taylor coefficients f(io,jo) of H(x,y) can be obtained via Cauchy's formulas. 

The starting point for proving Theorem 4 is the boundary value problem that satisfies 
the function H(x,0) (see Section 2.5). In Section 4.1, we show that H(x,0) is necessarily 
of the form S(w(x))/L(x,0), where S is a rational function. In Section 4.2, we prove that 
S(X) is actually a polynomial of the first order, and in Section 4.3 we find the expression 
of this polynomial (i.e., the constants \i and v in (4.2)). 

4.1. Expression of the solutions to the boundary value problem. For a given non- 
real zq and any non-zero harmonic function f(io,jo), consider the function 

{ ' L(x, 0)H(x, 0) -Zq 

We first notice that for xoaX([j/i,l]) (the boundary of Sfy), 7i(x) ^ oo. Indeed, remember 
(Section 2.5) that L(x,0)H(x,0) takes real values on X([yi,l]) and that zq is non-real. 
Further, let us denote by N p the number of poles that %{x) has within &x- First, by the 
principle of isolated singularities [24, Section 3.2], we obtain that N p < oo (for any non-real 
zq). Second: 

Lemma 5. For any non-real zq, we have N p 1. 

Proof. Assume that N p = 0. Then T~L(x) satisfies (i), (ii) and (iii), so that by Lemma 3, we 
conclude that T~L(x) is a constant function C. Accordingly, H(x,0) = {zq + 1/C)/L(x,0). 
Having zq + 1/C ^ is impossible, as 1/L(x, 0) is either not a power series (if L(0, 0) = 0), 
or is a power series whose certain coefficients are negative (if L(0,0) 7^ 0). It follows that 
zq + 1/C = 0, and thus that H(x,0) is identically zero. Accordingly, f(iQ,jo) = for any 
^0)j0 ^ 0, which is a contradiction with our starting assumption. □ 

Though this is not crucial for our analysis, it is interesting to note that we can prove, 
purely analytically, that N p is intrinsic, in the sense that it does not depend on zq: 

Lemma 6. N p is independent of the choice of zq G C \ R. 
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(4.3) — / = N Z -N P = -N p , 



Proof. The number of poles N p can be calculated as follows (see [24, Section 3.8.2]): if the 
contour X([yi,l]) is traversed in the counter-clockwise sense, and if N z (resp. N p ) is the 
number of zeros (resp. poles) of T-L(x) inside X([yi, 1]) (they are counted as many times as 
their order indicates), we have 5 

J_ f %'{x) 

2vri Jx([ yi ,i]) ti{x) 

since N z = (indeed, H(x,0) ^ oo for x £ &x)- By a classical reasoning (N p takes only 
integer values and is continuous w.r.t. zq in the positive (or negative) half plane), we obtain 
that N p is constant for zq in a given half plane. Since the coefficients of L(x, 0)H(x, 0) are 
real, we reach the conclusion that N p takes the same constant value for zq in the two half 
planes, and the proof of Lemma 6 is completed. □ 

We now show the existence of a function C(x) such that T~L{x) — C(x) satisfies (i), (ii) 
and (iii). With Lemma 3, it will then be obvious to deduce that T~L(x) — C(x) is a constant 
function. Let us call xi, . . . ,x q the poles of 7i(x) within &x, with respective multiplicities 
n\, . . . , n q , so that ni + • • • + n q = N p . Since the function 

1 

w(x) — w{x() 

has a pole of order 1 at xg (see property (W3)), is analytic in \ { x i} an d continuous in 
<Sx \ {xe}, we can find a polynomial Pi of order ri£ such that Pe(0) = and 

1 

w(x) — w{xi) 

is analytic at xg. The following choice for C(x) is thus suitable: 

c( x ) = m=i ^ ( K w( x )-w(x e )) ■ 

We conclude that the function T~L(x) — C(x) is constant. Since w(l) = oo (see (W4)), we 
have C(l) = (remember that P^(0) = 0). Furthermore, we have H(1,0) = oo (Lemma 
2), so that H(l) = 0, and the constant function %{x) — C(x) is zero. Finally, we have 

#M) = TT^V ( z + 



%{x) - P e 



L(x,0) \ u C(x) 

4.2. Form of the rational function. To summarize Section 4.1, given some zq G C \R, 
we have constructed a rational function 



such that H(x,0) = (zq + 1/R(w(x)))/L(x,0). We now prove that q = 1 in (4.4), and that 
the polynomial Pi has order 1. 

Lemma 7. The rational function 1/R(X) must be a polynomial. 



priori, this Cauchy's theorem cannot be applied as such, as the contour X([yi, 1]) is not included in 
a domain where H(x) is analytic. However, since H(x) can be analytically continued through any point of 
X([yi, 1]) \ {1} (see [21, Section 4] for a similar statement), and since H(x) — » oo as x — > 1, this result still 
holds. 
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Before proving Lemma 7, note its following consequence: 

Corollary 8. We have q = 1 in (4.4) and P\{X) = fJ,X ni , for some fi £ C and n\ € Z + . 

Proof. It is elementary to prove that for a rational function R(X) as in (4.4), a necessary 
and sufficient condition for 1/R{X) to be a polynomial is that q = 1 and Pi(X) = \xX ni 
in (4.4). □ 

Proof of Lemma 7. Assume that 1/R(X) is not a polynomial. Then there exists z\ G C 
such that \/R{zi) = oo. Introduce t\ £ Sfy such that 6 w(ti) = £i- Note that t± ^ 1, as 
w(l) = oo. We conclude that H(t\,0) = oo, which is an apparent contradiction with (i) 
or (ii) of Section 2.5. □ 

So far, we have shown that 

zq + n(w(x) - w(xi)) ni 



(45) H ^ L(x.O) 

Lemma 9. We have ni = 1 and (i£R \ {0} in (4.5). 

Proof. First, we have fj, ^ (we already saw that 1/L(x,0) cannot be the generating 
function of positive numbers). As x G (0, 1) goes to 1, both w{x) and H(x,0) are real and 
go to oo. For the polynomial zq + fi(X — w{x\)) ni to be real, we must have fiGR. If 
ni ^ 2, the same argument (together with the binomial formula) gives that w{x\) and Zq 
have to be real too. But zq € C \ R, so that m = 1. □ 

Finally, we have proved (4.2), and it remains to compute the constants fi and v. 

4.3. Computation of the constants. This part aims at computing the values of \x and 
v in (4.2) given in (4.1). 

Case L(0, 0) = 0. This is a simple case: taking the limit of (4.2) as x — > gives that v = 0. 
The value of \x then follows by examining the first term in the expansion of (4.2) at x = 0. 

Case L(0, 0) 7^ 0. A similar reasoning as in Section 4.1 and Section 4.2 gives that 

(4,) "(°->-i£±f 

where w(y) is the analogue of w(x), replacing the domain by £fy, the domain bounded 
by the curve Y"([o;i, 1]) and containing y\. Similar properties as (W1)-(W5) hold for w{y). 
Further, we fix w(y) so as ^(O) = and (3.3) holds (with exactly the same value of c). 

We show that for (4.2) and (4.6) to hold, we must have fj, = Jl. To this aim, evaluate 
the functional equation (2.4) at (Xo(y),y); we obtain 

= L(X (y), 0)H(X (y),0) + L(0, y)H(0, y) - L(0, 0)H(0, 0). 

Using (4.2) and (4.6), we may write 

(4.7) = n(w{X (y)) + + Jl(w(y) L(0, 0)H(0, 0). 



6 If ti belongs to the (open) set &x, then ti is unique, while if t\ is on the boundary, two values (namely, 
t\ and ti) are possible. 
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In the neighborhood of 1, we have Xq(u) = 1 + exp(i9)(y — 1) +o(y — 1), see [12, Equation 
(2.8)], so that 

, , c + o(l) c+o(l) 

= <i -*,<»))■>• = -(T^7»- 

Multiplying Equation (4.7) by (1 — j/) 7r ' 6 and letting y — >• 1, we obtain that /i = ju. 

Now we come back to (4.2) and (4.6). Evaluating them at x = and y = 0, respectively, 
we obtain that v = v = pi if(l, 1)/V- To find /i in the case L(0, 0) 7^ 0, we can use (4.7): 

pi,i/(U) 



The function in the denominator above is constant, and can be evaluated at any value of 
y, for example 0. 

4.4. Example. We consider the simple random walk, for which &x 1S the unit disc, see 
Section 2.4. With Theorem 4, in order to compute H(x,0), it suffices to calculate w(x), 
L(x,0), \x and v. First, the function 

= n X ^2 
(1 — x) z 

manifestly satisfies (W1)-(W5) (remember that 6 = 7r/2); see also (3.5). Furthermore, we 
have L(x,0) = x/4, so that L(0, 0) = 0, and for this reason v = 0. Finally, with (4.1), 
/i = /(l, l)/4, and thus 

This does match with (2.5). 



5. Applications 

In this section we present four applications of our results, which successively concern 
Martin boundary, (non-)algebraicity of the harmonic function, exit time from cones, and 
lattice path enumeration. 

5.1. Martin boundary. 

Brief account on Martin boundary theory [7]. For a (transient) Markov chain with state 
space E, the Martin compactification of E is the smallest compactification E of E for which 
the Martin kernels y 1— >• ky = Gy/Gy° extend continuously (Gy is the Green function, and 
we denote by xq any reference state). The set E\E is called the full Martin boundary. For 
a € E, x i—T- k% is superharmonic, and d m E = {a E E \ E : x 1— > is minimal 7 harmonic} 
is the minimal Martin boundary. Then, any superharmonic (resp. harmonic) function h 
can be written as h(x) = J-^kyfi(dy) (resp. h{x) = Jg E kyfi(dy)), where ji is some finite 
measure, uniquely characterized in the second case above. 



A harmonic function h is minimal if ^ g ^ h with g harmonic implies g = ch for some constant c. 
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Our contribution. In Section 4, provided that the drift of the walk is zero, we have proved 
the uniqueness of the harmonic function. In other words: 

Theorem 10. For all random walks satisfying (HI), (H2) and (H3), the Martin boundary 
is reduced to one point. 

In particular, all results [1, 2, 20, 22] regarding the uniqueness of the harmonic functions 
in some particular cases of walks in the quarter plane also follow from Theorem 10. 

5.2. Nature of the generating function H(x, y) and group of the walk. 

Nature of the generating function H(x, y). Let 9 be defined in (2.11). We have the following 
link between the nature (rational, algebraic, non- algebraic) of H(x,y) and 9: 

Proposition 11. The generating function H (x , y) is rational (resp. algebraic non-rational, 
non-algebraic) if and only ifir/6 belongs to Z (resp. Q, R\ Q). 

Proof. From the functional equation (2.4) and by symmetry, it suffices to prove the above 
statement for H(x,0). In fact, thanks to Theorem 4, it is enough to prove it for w(x). If 
tt/9 belongs to Z or Q, Proposition 11 easily follows from (3.4). If tt/9 6 R\Q, then it is 
a corollary of (3.3) that w(x) cannot be an algebraic function. □ 

A direct consequence of Proposition 11 is the following: 

Corollary 12. The harmonic function f(io,jo) is polynomial in %q and jo if and only if 



However, it is an open problem to decide whether f(io,jo) is homogeneous or not. For 
more details about this homogeneity problem, we refer to [22, Section 1.3]. 

(Ir) rationality of tt/9 and the group of the walk. As suggested by Proposition 11, tt/9 is 
indeed a key-quantity. However, it is not simple to check whether tt/9 is rational or not, 
and this is why we propose the following alternative way to find the (ir) rationality of tt/9: 
define Q(x,y) = x 2 y 2 L(l/x, 1/y), with L as in (2.3), and introduce 



It is quite easy to see that if Q(x, y) = 0, then Q(£,(x, y)) = Q(n(x, y)) = 0. In particular, £ 
and rj can be interpreted as automorphisms of the surface {(x, y) G CU{oo} : Q(x, y) = 0}. 
The group of the walk [11, Section 2.4] is then precisely the automorphism group W = (£, rj) 
generated by £ and 77. 

Further, it is proved in [12] that under (HI), (H2) and (H3), the group is finite if and 
only if tt/9 G Q. In other words, we can reformulate Proposition 11 by replacing u tt/9 £ Q" 
(resp. u tt/9 G R\ Q") by "the group W is finite" (resp. "the group W is infinite"), knowing 
that in many examples [3, 11, 18, 22], the order of W is already computed. 

polynomial P(io,jo) is homogeneous if there exists a £ Z+ such that P(\io, Xjo) = X a P(io, jo) for 
any iq, jo and A. 



tt/9 G Z. 
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Example. For the simple random walk, we saw that ir/9 = 2; therefore, the group W must 
be finite. Let us verify this fact directly: we have £(x,y) = (x, 1/y) and rj(x,y) = (l/x,y), 
so that W has order 4 (besides £ and rj, W includes the identity element, as well as the 
compound (oij = 77 o £ = (l/x,l/y)). Note that the order of the group is twice 7r/0, 
which is not surprising: in general, the order of W is 2inf{p € Z + \ {0} : p8/ir G Z}. 

5.3. Exit time from the cone. As written in the introduction, it is proved in [6] that for 
the zero-mean random walks in the quarter plane, Equation (1.1) holds, where V(io, jo) is 
some non-explicit harmonic function. Let f(io,jo) be the unique harmonic function found 
in Section 4. 

Corollary 13. There exists x > such that for any (io,jo) E Z5_, V(iq, jo) = x/(*0) Jo)- 

Remark. The conformal mapping tu(x) is related to the first exit time from the cone, as 
V(io,jo) governs the asymptotic tail distribution of the latter (Corollary 13). Actually, the 
link seems to be stronger, as the following holds: let F(n) and 6 be introduced in (1.1) and 
(2.11), respectively. It is proved in [6, Theorem 1 and Example 2] that F(n) = n~ n ^ 29 \ 
and in (3.3) that w(x) = (c + o(l))/(l — x) n ^ , with the same angle 0. 

5.4. Lattice path enumeration. A point in the lattice Z+ has eight neighbors, diagonal 
neighbors included. Choosing a subset of the set of neighbors and considering it as the set 
of admissible steps of a walk in the lattice, we obtain one of 2 8 possible models of walks 
in Zj_. Bousquet-Melou and Mishna [3] showed that, after eliminating trivial cases, and 
also those which can be reduced to the walks in a half plane, there remain 79 inherently 
different models. Recently, in [6, Section 1.5], the asymptotics of the number N n (x,y) of 
excursions (the number of walks starting at x and ending at y after n steps) was obtained: 
as n — > 00, 

N n (x, y) = CV{x)V'(y)p n n a {\ + o(l)), 

where the constants C, p and a are found explicitly, while V{x) and V'{y) are not. In 
fact, with the notations of Section 2, V(x) is the harmonic function for the random walk 
(X,Y), and V'(y) is the harmonic function for the reversed process —(X,Y). If the drift 
of the walk (see (H3)) is zero, then the drift of the reversed walk is zero too, and the 
results of Section 4 provide expressions for V(x) and V'{y). If the drift is non-zero, the 
harmonic functions have been found in [18]. They could also be obtained from Section 6 
below, where we present some extensions (for instance, to the non-zero drift case) of our 
approach and results. 

6. Extensions 

Throughout this paper, we mostly assumed that the random walks had zero drift (H3). 
However, our approach can also be applied to obtain the Martin boundary in the non-zero 
drift case (Section 6.1), and to compute t-harmonic functions (Section 6.2); we also say a 
word on a possible generalization to random walks with more general jumps (Section 6.3). 
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6.1. Random walks with non-zero drift. In this section, we consider random walks 
satisfying (HI) and (H2) (but not (H3)). We first recall the existing results in the literature 
on Martin boundary [16, 18]. Then we show how the approach in this paper allows us to 
retrieve these results. Finally we present an interesting phenomenon on the convergence 
of the harmonic functions of the non-zero drift case to the unique harmonic function in the 
zero drift case. In this section, we focus on ideas, and therefore just sketch the proofs. We 
illustrate all three paragraphs by the example of the simple random walk. 

Existing results. For random walks with non-zero drift killed at the boundary of the quarter 
plane, the minimal Martin boundary (see Section 5 for a brief account on Martin boundary 
theory) is homeomorphic to [0, 7r/2], see [16, 18]. In concrete terms, this means that the 
set of all minimal harmonic functions is homeomorphic to [0,7r/2], and therefore can be 
parametrized by 7 G [0, 7r/2]. Further, an expression for the harmonic functions / 7 is given 
in [18], in terms of solutions to certain boundary value problems. While this expression is 
typically complex, it may admit nice simplifications (in fact, when a certain group related 
to the walk is finite [18]). 

For instance, for the simple random walk, we have (see [18, Section 5.1]) 

r ( s+ (o) io - s-(oy°)j s+(o) jo if 7 = 0, 

(6.1) Mio,jo) = { (M7y°-M7) io )(M7y°-S-(7) i0 ) if 7G(0,tt/2), 

[ i s+(vr/2) Jo (? + (7r/2p° -?_(vr/2p ) if 7 = vr/2, 

where 



r( 7 ) 
r(7) 



1 ~ ~ (1 ~ tan 2 7)(l - 4p ,-iPq,i + 4p_ li op li o tan 2 7) 

2pi )0 (l - tan 2 7) 
1 ~ \/l - (1 - tan~ 2 7)(l - 4p„i ; oy>i i0 + 4p _iPo,i tan" 2 7) 
2pi n(l - tan" 2 7) 



and 



s± ( 7 ) = r ( 7 ) ± y/r(7) 2 -p_i )0 /pi i0 , s±(t) = r {l) ± \/ r {l) 2 ~ Po ,-i/Po,i- 



Retrieving the existing results. In this paragraph, we still consider non-zero drift random 
walks, and we explain how our methods provide the Martin boundary and even expressions 
for harmonic functions. In the zero drift case, the point 1 plays a crucial role, in particular 
for the boundary value problem of Section 2.5, see (i). This is a consequence of the radius 
of convergence (equal to 1) of H(x, 0) — or equivalently of the non-exponential growth of 
f(io, 1) as iq — > 00. For non-zero drift walks, the growth of f(io, 1) is now exponential. 
In fact, we could prove that this exponential factor necessarily takes the form 1/p, with 
p 6 [x2, X(y.2)] 9 , see (6.1) for the case of the simple random walk. We thus have infinitely 
many boundary value problems, parametrized by p, as follows: 

(i) H(x,0) is analytic in Wx, except at p E [x2,X(y2)] } where it has a pole 10 ; 

(ii) H(x,0) is continuous on &x \ {p}; 



^In the zero drift case, we have X2 — yi — X{y2) = 1. In the non-zero drift case, the points X2 and j/2 
are defined as in Section 2.4, by replacing the inequality (2.9) by the following: |xi| ^ X2 ^ 23 ^ \ x i\ an d 
|yi I ^ V2 < Wa < \va\- 

In fact, the pole is simple if p € [X2, -X^jfe)), and has order two if p = X(y2). 
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(iii) For all x on the boundary of L(x,0)H(x,0) - L(x,0)H(x,0) = 0. 

The resolution of (i), (ii) and (iii) (that we shall not tackle here) would use an analogue of 
the conformal mapping w, namely, a function such that: 

(Wl) w is analytic in except at p E [x2, X(y2)], where it has a pole 10 ; 
(W2) w is continuous on &x \ {p}', 

(W3) w establishes a conformal mapping of <£x onto the complex plane cut C along a 
segment; 

(W4) For all x on the boundary of w(x) = w(x). 

At this point, we wish to emphasize the following main difference between the zero drift case 
and the non-zero drift case: the expression of w(x) becomes more elaborate (more precisely, 
in place of sine and arcsine functions, it now involves two p-Weierstrass functions with 
different periods 11 ). Despite this intrinsic complexity, expressions for w(x) are available 
in the literature: see [18, Section 3.1] for p = X2, and [21, Section 2] for simplifications of 
w{x) in some particular cases. Once w{x) is known, Equation (4.2) holds, and an explicit 
expression for H (x, 0) follows. 

In some sense, this approach is more efficient than the one used in [18], since we give 
here an expression for the harmonic functions directly in terms of the function w(x), and 
not in terms of solutions to some other boundary value problems. 

For the simple random walk, <£x is the circle with center at and radius \/pi,o/P-i,o- 
The function 

W X (x -p){x- (pi,o/P-i,o)/p) 
therefore satisfies (W1)-(W4), and with (4.2) we conclude that 

w(x) l/po,i (l/p ,i)(P-i,oM,o) 



H(x,0) 



L(x,0) (x -p)(x - (pi,o/p-i,o)/p) (1 - s+(7)z)(l - s-irfx) ' 

provided that 7 is defined by l/s+(7) = p. The latter equation has a solution if and only 
if p 6 [x2,X(jj2) = y/pifi/p-i,o}- Then, using the functional equation (2.4), this forces the 
pole of H(0,y) in Sfy to be Yq(p), that is, l/s + (j). 

Convergence of the harmonic functions of the non-zero drift case to the harmonic functions 
in the zero-drift case. Consider a random walk without drift, and with jumps {pij}— i<i i^i- 
Denote by / its unique harmonic function. Let us introduce any family of random walks 
with non-zero drift and jumps {^"j}— such that for all i and j, linin-^oo q^j = Pi,j- 
For any n, there exist infinitely many harmonic functions parametrized by 7 £ [0, vr/2]. 
We normalize all these functions by requiring that /(l, 1) = /"(l, 1) = 1. The main result 
in this paragraph is the following convergence result: for any io,jo and any 7, we have 

(6.2) lim / 7 n (i ,io) = /(*o,jo). 

n— too 1 

We now mention two open problems concerning (6.2): 



11 This is due to the fact that a certain Riemann surface — namely, {(x, y) G (C U {oo}) 2 : L(x, y) — 0}- 
sees its genus passing from to 1. 
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• The proof that we could write for (6.2) is purely analytical. It is an open problem 
to decide whether this can be obtained by using only general harmonic analysis. 

• The convergence (6.2) to the same harmonic function f(io,jo) can be proved 
independently of knowing that in the zero drift case, there exists a unique harmonic 
function. It is an open problem to prove that the convergence (6.2) implies the 
uniqueness of the harmonic function in the zero drift case. 

Let us now illustrate (6.2) with the simple random walk, for 7 £ (0, 7r/2). With (6.1), 
the function H y (x,y) = Ei jo»l A^ ' ioK° _ V 0_1 equals 

(«+(7) ~ g-(7))(g+(7)^-S r -(7)) 

(1 - xs+(j))(l - xs-(i))(l - y?+(7))(l - ys-(j)) ' 

As the drift goes to 0, the four points s+(7), £-(7), s— + (7), £-(7) all converge to 1, so that 
dividing by the constant (s + ( r y) — S-(j))(s+( / y) — £-(7)), the generating function H^{x,y) 
goes to 1/[(1 — x) 2 (l — y) 2 ], which indeed is the generating function for the zero-drift case, 
see Section 2.3. The convergence of the generating functions yields the convergence of its 
coefficients. 

6.2. i-harmonic functions. In this work, we have determined the harmonic functions, 
that is, the i-harmonic functions (satisfying Pf = tf) for t = 1. It is also interesting to 
consider i-harmonic functions other values of t (many motivations are given in [26, Chapter 
2]). It turns out that our approach can be extended to find the latter: indeed, the key 
functional equation (2.4) still holds, if now 

L(x,y) = L(x,y;t) = xy[Y^- 1<id<1 Pi,jX~ l y~ J - t]. 

Accordingly, the methods presented in this article can be used to find t-harmonic functions, 
with, nevertheless, the following additional difficulty concerning the conformal mapping 
w(x): mutatis mutandis, as for the non-zero drift case (see Section 6.1), introducing the 
variable t would have the effect of involving p-Weierstrass functions instead of sine and 
arcsin functions in the expression of the crucial conformal mapping w(x). 

6.3. More general jumps. In this part we consider random walks in the quarter plane, 
that do not satisfy the small step hypothesis (HI). Let us only assume that it has bounded 
jumps, which means that there exist and J^, such that pij = if i > I + , i < I—, 
j > J + or j < J~ . We are interested in harmonic functions satisfying to (PI), (P2) 12 and 
(P3). Let us define the sectional generating functions 

F j0 (x) = E^i/OWo)** - 1 , O(y) = EjozifihM - 1 . 

The following functional equation then holds: 

(6.3) L(x,y)H(x,y) = £ lcjo<J+ F^(x)y^ 1 T.i-^i+Pi^-^y'^ 

~ Ei<i </+ /(*o.io)x i °-V' - 1 Ei ^j+ Pijx-^y-ti- 1 ). 

12 In fact, (P2) has to be replaced by f(io,jo) = if io ^ or jo ^ 0. 
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This equation is far more complicated than (2.4), because instead of 2 unknowns (namely, 
H(x, 0) and H (0, y)) in the right part of (2.4), there are now I + + J + unknowns (namely, 
F 30 {x) for 1 ^ jo ^ J + and G l °{y) for 1 ^ io ^ / + ) in the right-hand side of (6.3). To our 
knowledge, it is an open problem to solve (6.3) in general. 

However, if we restrict ourselves to I + = J + = 1 (with I~ and J~ arbitrarily large), 
then (6.3) is nothing else but (2.4), as F 1 (x) = H(x,0) and G 1 {y) = H(0,y). In particular, 
the approach developed in the present paper can be applied for any such random walk, 
provided that a conformal mapping w(x) can be found: if I~ or J~ are different from — 1, 
then the kernel L(x, y) becomes more complicated, and it is an open problem to construct 
a suitable conformal mapping w{x) with properties (W1)-(W5). 

Appendix A. Proof of Lemma 2 

By symmetry, it is enough to consider only H(x,0). We shall prove that the radius of 
convergence of H(x, 0) is at most 1. This turns out to be sufficient, since Lemma 3 directly 
implies that H(1,Q) = oo. To prove that the radius of H(x,0) is smaller than 1, we show 
that limsup i _ >00 (l/i) log h(i, 1) ^ 0. 

Let us define Ti 1 j 1 = inf{n ^ : (X(n),Y(n)) = the first hitting time of state 

(h,jl), as well as the Green functions G^-° = J2n>o P (io,io) i(. X ( n )> Y ( n )) = (hji)]- The 
classical Harnack inequality yields 

f(i,l) 1 _ Glj 

/(l,l)^P (1|1) [T i>1 <cx ) ]~Gy- 

First, we have that sup,^ G*'} < oo. Indeed, the Green functions G-'} for the quadrant 
are smaller than the Green functions G l A associated with the half plane Z x Z+. The 
latter domain being invariant by horizontal translations, the Green functions G\'\ do not 
depend on i. Accordingly, we have (see [13] for the last inequality) 

supG-'} ^ supG-'} = Gl'l < oo. 
To conclude, it is enough to prove that 
(A.4) liminf-logG^ 1 > 0. 

To that purpose, let us first use the last inequality in [13, Proof of Proposition 4.2] (which 
is valid not only for the half plane — the topic of [13] — but also for the quarter plane); we 
obtain 

(A.5) lim inf - log G X A ^ - inf L T] (0), 

«->oo I ' <j) 1 ' J 

where the infimum is taken over all trajectories 4> : [0, T] -)• such that 0(0) = (0, 0) 
and <p(T) = (1,0), and where the rate function i]o,T](0) equals L((j)(t))dt, with 

L(v) = {\ogR)*{v) = sup ((a,v) - log R(a)), R{a) = ^ p itj exp((a, (i, j)}). 
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Obviously, inequality (A. 5) still holds if one takes the infimum over all linear paths (ft, i.e., 
(ft(t) = (t/T)(l,0) for t € [0,T]. For such functions (ft, we have 

inf/ [0>r] (</))=T-(log J R)*(l/T,0) 

Then, as (see [23, Page 35]) 

mf(r-(lo gj R)*(l/T,0)) = 0, 
we reach the conclusion that (A. 4) holds. The proof of Lemma 2 is completed. 
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